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1. Introduction 

The Standard Model (SM) has been completed with the discovery of the Brout-Englert-Higgs 
(BEH) boson [p at the Earge Hadron Collider (EHC) 0. Besides, the outstanding performance 
of the machine and of the EHC experimental collaborations have also allowed to translate the non¬ 
observation of any significative departure from the SM predictions into stringent limits on New 
Physics (NP). What, in particular, implies that if NP is going to show up at future runs, it must 
predict relatively small cross-sections into final sfafes already scrutinized. Obviously, small cross- 
secfions can be always obfained banishing fhe NP scale fo higher energies or coupling if fo fhe SM 
weakly enough. The relevanf quesfion is if fhere is room leff for new inferacfions near fhe TeV scale 
wifh effecfive couplings of af leasf Elecfro-Weak (EW) sfrengfh. Whaf also brings fo fhe quesfion 
of whefher NP can manifesf af fhe Infernafional Einear Collider (lEC) if no deparfure from fhe SM 
can be esfablished af fhe EHC. 

Specific SM exfensions can predicf some new parficles fo be in pracfice invisible af fhe EHC 
even if fheir masses are near fhe EW scale. Indeed, fhe EHC discovery reach can be reduced 
allowing fhese heavy particles fo decay info new channels wifh larger backgrounds, making in any 
case fhem wider and evenfually undefecfable. Their producfion, and fhen fheir defection, can be 
also made more difficull if fhey are consfrained fo be pair-produced by imposing fhe corresponding 
selection rule. There are many of fhese examples in fhe liferafure. Eor insfance, sfealfh gluons [||] 
provide an example of fhe firsf case, and supersymmefric models wifh R-parify 0 of fhe second 
one. Offen fhese properfies reducing fhe EHC discovery pofenfial are also refained af fhe lEC, nof 
being fhen fhis machine an alfernafive for fhese NP searches. Alfhough fhe acfual sensifivify 
for a given process will depend on which circumsfance plays a more imporfanf role, ifs lower Cenfer 
of Mass Energy (CME) or ifs cleaner environmenf, and evenfually fhe available beam polarization 
af fhis lepfonic machine [^. 

One can be more drastic, however, and assume fhaf fhe new parficles do nof couple fo quarks 
and gluons, nor fo EW bosons, buf only fo SM lepfons, af leasf fo sfarf wifh. Trying in fhis way fo 
reduce as much as possible fhe EHC discovery reach, buf nof fhe lEC one. Nonefheless, even in fhis 
exfreme scenario fhe EHC seems fo be able of partially settling fhe question, whenever elecfrons 
(posifrons) are nof relevanf. A defailed sfudy of fhe EHC and fhe lEC pofenfial for fhe discovery 
of fhese lepfophilic inferacfions is presenfed in Ref. 0. We summarize fhe main resulfs presenfed 
fhere in Section buf concenfrafing on fhe parficular case of an exfra gauge boson coupling fo 
muon minus fau Eepfon Number (EN), |0], as example. We fhen compare fhe EHC and lEC 

pofenfials wifh fhe pofenfial of a Eufure Circular hadron Collider (ECC) wifh a CME ^/s = 100 
TeV [^] ^ in Section Alfhough such a lepfophilic gauge boson fends fo mix little wifh fhe SM 
Z boson, fhe kinefic mixing can be large in generic models [jlT]] and fhen, fhe Z^_^ confribufion fo 
fhe Drell-Yan process, qq ^ ^ ll, could be evenfually sizable af a large hadron collider. We 

discuss fhis scenario in Section ^ finding fhaf fhe limifs on fhe complemenfary process af fhe EEP 

2, —>• Z^_.j. —)■ U, already consfrain fhe kinetic mixing fo be small for relafively low vector 
boson masses of few hundreds of GeV. However, Drell-Yan producfion af fhe EHC (ECC) provides 
fhe mosf sfringenf bound on such a mixing, especially for larger masses of fhe order of fhe TeV. 

'Future Circular lepton Colliders have also deserved consideration |^], as well as a Large Hadron electron Collider 
(LHeC) [p^] with different options. 
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Figure 1: Leading diagram for Z^_t. production at the LHC, for negligible kinetic mixing. 

Thus, four-lepton production in Fig. |I] is the only leptophilic final state available for vanishing 
kinetic mixing but qq —>• 7 ,Z,Z^_^ —)• /fallows to set quite strong limits on these interactions for 
large kinetic mixing. Final comments are collected in the Conclusions. 


2. Limits on at the LHC and the ILC 


As we have recently argued in Ref. 0, the only leptophilic particle coupling to SM lepton 
pairs in a renormalizable way is a new neutral vector boson Z'. If this boson was a scalar, it 
would have to be the component of a multiplet transforming as the BEH one or including a doubly- 
charged component, too (see Ref. |[T^ for a detailed discussion of the lowest order couplings and 
production mechanisms for such a scalar). Thus, implying in either case that it would couple to 
SM gauge bosons with EW strength and hence, it would not be a leptophilic particle. (Although 
at the end, leptophilic additions are expected to partially lose their character because higher order 
effects may result in a small mixing between the different sectors of the theory, thus coupling to 
EW gauge bosons and to quarks, too, but with suppressed couplings.) Bosons with spin higher than 
1 do not have renormalizable couplings to lepton pairs, either. The generic Eagrangian describing 
the corresponding vector boson interactions reads: 


.jZ’z' = - , (2.1) 


where g'^^ are arbitrary dimensionless couplings to the SM Eeft-Handed (EH) and Right-Handed 


(RH) lepton multiplets, L,l = and /*, respectively, with / = labeling the charged- 

lepton family. Even though we allow for arbitrary flavor and chiral interactions, the EW gauge 
symmetry remains unbroken because the new couplings of the EH charged leptons and their neu¬ 
trino counterparts are equal. However, the Yukawa couplings giving masses to the charged leptons 
in the SM do not need to preserve such a hypothetical gauge symmetry. As a matter of fact, if a lep- 
tonic Yukawa term has to be invariant under the new symmetry, the scalar field coupling to the pair 
of EH and RH leptons must have a commensurate charge adding to zero with the lepton charges. 
Thus, in general, the existence of a new leptophilic interaction described by Eq. ( ^T| ) requires 
non-minimal scenarios with a non-standard lepton mass generation and/or extra scalar multiplets 
for the model to be realistic. 

An example of the first scenaiio, without extra light scalar multiplets and hence, without 
Yukawa terms involving leptons (because the SM BEH doublet must be a singlet under the new 
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symmetry in order to quarks get SM masses), is the composite Higgs model with tau custodians in 
Ref. [jlsl]. In this case, once the extra heavy vector bosons present in the model acquire a mass, 
lepton Yukawa couplings are generated at higher orders in perturbation theory. 

A prime example of a model with an extra leptophilic interaction and possibly extra scalar 
doublets is the combination of muon minus tau LN, — L^, ^ with diagonal couplings 
given by the charges 


Multiplet 

‘-(r) 

Mr 

‘“4*) 

Fr 

(2.2) 

Charge 

1 

1 

-1 

-1 



times the gauge coupling constant g'. This is the only gauge addition to the SM with vanishing 
couplings to electrons and anomaly free in the absence of extra fermions li- On the other hand, this 
model also allows for a good estimate (which can be also easily extrapolated [|^) of the LHC (and 
the ILC) potential for the discovery of leptophilic interactions. Its reach appears to be relatively 
reduced (in either case) because the two-lepton signal of the new resonating vector boson must in 
principle emerge from a four-lepton sample. This means that we must deal with relatively small 
cross-sections of few fb at the LHC for = 14 TeV (and similarly at the ILC) and hence, that there 
is only enough sensitivity for an eventual discovery when the new gauge couplings are relatively 
large and/or the Z' masses are somewhat small. This then excludes possible electron signals if the 
limits from LEP 2 must be fulfilled I®- Assuming universality, for these limits do not show a 
large flavor dependence, the corresponding 95 % Confidence Level (C.L.) bounds read |]_5] (see 
also [[T^): 


tee 

— < 0.12,0.16 TeV^i 
Mr 


(2.3) 


being too stringent for observing a leptophilic vector boson in a four-lepton final state at the LHC 
(ILC). We assume that the model gauging L^ — L^- also includes a SM singlet scalar giving a 
relatively large mass to the new gauge boson. The originally proposed realization addressed the 
possibility of a relatively light Z'^_^. In our case, the extra vector boson can be made arbitrarily 
heavy. We will advocate for the extra scalar doublets transforming non-trivially under the new 
gauge symmetry only to generate an appreciably mass mixing with the SM gauge bosons, when 
needed. Neutrino masses are very tiny and we assumed they are generated at higher orders in a 
non-minimal scenario. In any case, in the following the parameters used to fix the model in the 
Monte Carlo simulations shall be only the new gauge coupling g' and the heavy vector boson mass 
M 71 ; while the total width writes: 


Tyt = 


g'^ 

—Mr 

4k 


(2.4) 


where we neglect any Z'^_^ decay into new scalar pairs. 

The LHC and ILC limits on a generic leptophilic vector boson only coupling to muons and taus 
are thr'oughly discussed in Ref. [^. Let us review here the case of the Z'^_^ above. In general, such 

^In order to give a large enough mass to the new gauge boson, the model must also include at least one extra scalar 
field transforming trivially under the SM gauge symmetry but not under — L^. 
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Figure 2: Discovery and exclusion limits for as a function of its mass at the LHC with ^ = 14 TeV. 
On the left (right) we plot the bounds obtained from 3/4 plus missing energy (4/i) events. We also draw the 
bounds from neutrino trident production (straight line), for comparison (see the text for details). 

a vector boson will show up in the four-lepton final states. As a matter of fact, the best limits at the 
LHC can be generally derived from the 3/4 plus missing energy events. Although for this model the 
bounds from 4/4 are somewhat similar. In Fig. ^ we plot the corresponding discovery and exclusion 
bounds for the 3/4 plus missing energy (left) and 4/4 (right) samples. As explained in Ref. [|6|], these 
curves can be obtained from the plots in this reference multiplying by the appropriate factors. In 
particular, the exclusion bound can be read from Fig. 3 with ^ 1 there multiplying 

the corresponding coupling constant for a given Z' mass by s/l, to take into account the effect of 
the extra vector boson decay channel into taus. The event generation (codes used) and selection 
(applied cuts) are described in that reference, too. For comparison, we also show the bound from 
the neutrino trident production (straight line) [|T7|], taking into account CHARM-II and CCFR data. 

Whereas only signals involving muons can emerge at the LHC for the cross-sections at 
hand, as pointed out in Ref. [^, the ILC could observe the mixed channels with tau leptons. Thus, 
in Fig. ^ we show the corresponding discovery and exclusion limits but for 4/4 and 2/42 t events, 
respectively. Again, the bounds can be read from that reference. In particular, the limits obtained 
from 4/4 events can be read from Fig. 7 (right) there multiplying the corresponding coupling 
constant for a given Z' mass by s/l, to account for the extra decay channel into taus. Details on the 
event generation and selection cuts are given in the same reference. 

As emphasized in Ref. [|6|], the LHC can be sensitive to a with a mass up to ~ 1 TeV for 
g' of order 1; while the ILC is sensitive to Z^_^ masses below 300 GeV but it can eventually allow 
for a direct measurement of the couplings to tau leptons. 

3. Limits on at a future hadron collider with ^/s = 100 TeV 

As already pointed out, the LHC and the ILC will improve the limits on a new leptophilic Z' 
unmixed with the SM, but for relatively large coupling constants and relatively small vector boson 
masses. Their reach will be appreciably enlarged at a larger hadron collider with a CME of y/s = 
100 TeV, however. Similarly as for the LHC in the previous section, in Fig. ^we plot the discovery 
and exclusion bounds for the 3/4 plus missing energy (left) and 4/4 (right) samples. The codes used 
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Figure 3: Discovery and exclusion limits for as a function of its mass at the ILC with ^ = 500 GeV. 
On the left (right) we plot the bounds obtained from 4)i (2ji2z) events. The constraint from neutrino trident 
production (straight line) is also shown for comparison. 


channel 



4/i channel 



Figure 4: Discovery and exclusion limits for as a function of its mass at the future hadron collider 
(FCC) with ^/s = 100 TeV (thick curves). On the left (right) we plot the bounds obtained from 3/4 plus 
missing energy (4/4) events. We also draw present bounds from neutrino trident production (straight line) 
and the LHC limits for y/s = 14 TeV (thin curves), for comparison. 


for the event generation are the same as in Ref. |j^. The selection cuts applied to the 3/4 plus 
missing energy sample must be optimized, however, to obtain similar limits to those derived from 
the 4/4 channel. This is necessary to get rid of the background events in that case, too. Thus, we 
increase the basic cut on the muon transverse momentum in Table 1 in the previous reference from 
50 to 100 GeV; and the cuts on the missing transverse energy and on the transverse mass from 100 
to 150 GeV and from 110 to 140 GeV, respectively. As a result, the FCC bounds on the gauge 
coupling improve by almost a factor of 5 those at the LHC for low masses; whereas its mass 
reach is near 2.5 TeV for g' ~ 1. 


4. kinetic mixing with the SM and di-lepton versus four-lepton signals of a 
leptophilic vector hoson 

There is a wide literature on the implications of the kinetic mixing between abelian gauge 
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factors [18] in SM gauge extensions [p^]. We will follow Ref. [ pi] ] because it emphasizes our main 
observation in this section concerning the size of the kinetic mixing of an extra abelian interaction 
with the SM: the effective mass mixing between a new Z' and the SM Z boson must be rather 
small, as required by present data, but the kinetic mixing can be relatively large because it only 
enters alone in the Z' interactions. Indeed, all deviations from the SM predictions for processes 
involving only SM particles and, in particular, testing the Z boson interactions are proportional to 
the former and hence, it must be tiny (below the per mille level) because no significant departure 
from the SM has been observed. While a non-vanishing kinetic mixing modifies the Z' current 
adding an extra term to it, linear in the mixing and proportional to the SM hypercharge, allowing to 
constrain its size in an independent way in processes exchanging the new vector boson. In the case 
of the Z'^_^ it means, for example, that although this gauge boson is leptophilic to start with, it also 
couples to quarks because they have non-zero hypercharge but with couplings suppressed by the 
kinetic mixing. The relevant question is then how small this mixing must be on phenomenological 
grounds. 

Let us be quantitative. Following (the notation of) Ref. o the effective mass mixing, 
rotating the current eigenstate gauge bosons to the heavy mass eigenstate basis, can be written 
(note that we used ^ before to specify the RH Z' couplings in Ref. [0]): 


tan2i§ Ri —2 


5M^ + sin;^ syv M| 

M, 


(4.1) 


where 8 M^ is the mass term mixing the gauge boson current eigenstates and sin;^ the corresponding 
kinetic mixing (and sw the sine of the EW mixing angle). 


^.ixing ~ 5M^ Z'^_, Z^ - ^ z;% . (4.2) 

We only keep the leading terms in the limit of ^ and % small (fields are in the physical basis). We 
will concentrate on the case <C Z> besides, which can be the result of a large mass hierarchy 
M| <C , as we have formally required when deriving Eq. or of a partial cancellation of 
the two terms in the numerator of this equation, or of both. A large Z'^_^ mass can be generated by 
a large vacuum expectation value (vev) of a scalar field transforming trivially under the SM gauge 
symmetry; whereas the mass mixing cancelling the kinetic term can be obtained from a scalar 
multiple! transforming non-trivially under both the SM and the extra abelian gauge symmetry and 
developing a vev of the appropriate size. ^ 

The new gauge boson interactions with the SM fermions are described by the Eagrangian 

^z' , « - sinz — Jl) 4% 

= -[(g' + sinZ;^)(v^L7cTV^L + AtL7cTAtL) + U' + sin^ —)AtR7cTAtR 

Z Cw Cw 

^Many phenomenological analyses neglect kinetic mixing effects, as, for instance, when determining the ZZ' mixing 
at the LHC [ po| | ; or the mass mixing generated by the vev of scalar fields transforming non-trivially under both abelian 
factors, as in Higgs-portal studies [0] . Although the kinetic mixing can be large and the effective mass mixing small in 
specific models [p^. 


7 








Leptophilic interactions at large colliders 


Francisco del Aguila 




Figure 5: On the left we plot the total cross-section for pp —^ for different values of sin;if (see the 

text) as a function of the new gauge boson mass at the LHC and at the FCC. On the right we show in the 
g' — M^i plane the excluded (upper) regions by the LEP 2 and ATLAS bounds (see the text). 
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+(-g + smX:r^){VrLYcjV-ch + Tl/ctTl) + {-g' + sin;i:—) tr7ctTr 

+ sin X ( v^l Ja v^l + ) + sin ;i: —Fr Jae^ 

^ _ 2 ^ c _ 

-sin;i^--(ML 7 CTnL + t/L 7 CT^^L) -sin;|^--nR 7 (TnR + sin;|^-- dRyadK +.. 


(4.3) 


6 c^/ 


3 C]v 


3 civ 


■]Z[ 


p-T 


where the dots stand for the same quark couplings as for the first family but for the other two quark 
generations. Thus, a non-vanishing kinetic mixing modifies fhe currenf, adding an exfra 
confribufion linear in sin;^ and proporfional fo fhe SM hypercharge, Y = Q — T^, wifh Q fhe elecfric 
charge and fhe fhird componenf of isospin. Therefore, if adds couplings fo elecfrons and fo 
quarks, implying a new confribufion fo Drell-Yan producfion af fhe LEP 2, —)• Z^_^ —)■ U, and 

af fhe LHC, qq —)• Z^_^ —)■ ll. In fhaf case, presenf limifs on fhe former also resfricf fhe allowed size 
of fhe kinefic mixing because no deparfure from fhe SM predictions has been observed buf only for 
relatively low Z'^_^ masses. In facf, fhe limifs in Eq. {13) franslafe info (using fhe couplings in Eq. 
( 0 ): 


< 0.07 TeV-2 , 

Ml ^ 


M, 


< 0.45 TeV^ 


z' 

^ii-x 


ieeee) . 


(4.4) 


These bounds apply fo fhe RH operators, which provide fhe mosf sfringenf limifs due fo fhe chi- 
ralify dependence of fhe hypercharge couplings. The lepfonic couplings in Eq. ( |4.3| ) are also non- 


universal, buf lepfon universalify was assumed in order to derive Eq. ( |2.3| ). In fhis sense, allhough 
if is a ralher good approximalion, fhe quanlilalive analysis is approximate. We also neglecf higher 
order ferms in sin;^. In any case, allhough Ihese bounds resull in small gauge couplings to quarks 
for small gauge boson masses, fhe direcl Z'^_^ production cross-section can be slill sizable af fhe 
LHC if fhe LEP 2 limifs on sin x are saluraled. As a mailer of facl, Ihese bounds become unreslric- 
live for large veclor boson masses. To be quanlilalive, we plol in Pig. (lefl) fhe folal new vector 
boson production cross-section a {pp —)• Z^_^), which is proportional to | sin;^|^, af fhe LHC and 
fhe PCC for sin;^ = 1 (dashed lines) and for sin;^ saluraling fhe firsl inequalify in Eq. ( |4.4[ ) wifh 
g' = 0.25 (solid lines). (Obviously, we always require |sin;^| < 1.) The second inequalify, which 
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is independent of g', is less restrictive in this case, and the corresponding curve lies in between the 
other two lines plotted. As it is apparent, the exclusion of two-lepton cross-sections below the fb 
by the LHC and the FCC will provide quite stringent limits on the kinetic mixing as well as on the 
mass and coupling of a As a matter of fact, the limit reported by the ATLAS Collaboration 

at the LHC with ^/s = 8 TeV on the production cross-section of a new neutral gauge boson decay¬ 
ing into two muons a {pp —)• Z') Br (Z' —is already more restrictive, for a generic 
Z^_^ — p^p^ branching ratio of 1/3, than the LEP 2 one, as shown in the same picture (lowest 
curve). It must be again emphasized that this figure is only for illustration, because we assumed a 


small enough kinetic mixing when deriving the Lagrangian in Eq. (4.3), for instance, but draw the 
curves for arbitrary sin;^. Similarly, the gauge boson branching ratio into muons may be smaller 
than 1/3 for non-zero kinetic mixing. In fact, the ATLAS bound defines a band depending on fhis 
branching rafio (neglecting fermion masses): 


Br (Z' 




^p^p ) = 


1 -|-1.5 V-1-0.625 


3 + 5x2 


X = 


sin/ e 
g' civ 


(4.5) 


which ranges from 0.01 fo 0.45. Wifh fhe same warning, on fhe righf panel we show fhe excluded 
(upper) regions in fhe g' — plane by fhe firsl LEP 2 inequalify in Eq. ( [4.4[ ) for sin x safuraf- 
ing fhe second inequality in fhe same equafion (lowesf dashed line) and for fwo ofher illusfrafive 
values of sin;^ = 0.1,0.01. The smaller fhe kinetic mixing is, fhe larger fhe allowed region in fhe 
g' — Mz'^ ^ plane is. We also draw fhe ATEAS Collaboration limif (solid line) on muon-pair pro¬ 
duction mediated by a new neufral gauge boson [|^] safurafing fhe firsl inequalify of Eq. (Q, for 
illusfralion. Obviously, fhe observation of a Z'^_^ resonance in four and fwo-muon evenfs would 
allow for a determination of bofh paramefers, g' and sin;^, once fixed fhe gauge boson mass. 

The previous discussion is based on phenomenological grounds and hence, a final comment on 
the expected size of the kinetic mixing is in order. As the studied model is renormalizable, quantum 
corrections can be actually calculated, being in general the renormalization of the kinetic mixing 
quite small and below the per mille level. Thus, a relatively large mixing as the one saturating 
Eq. ( |4.4D requires a similar value at the matching scale. Therefore, although Drell-Yan production 
overcomes four-lepton production except for rather small kinetic mixing, the latter provides the 
most stringent model-independent limits for gauge boson masses below a few TeV. 


5. Conclusions 

We have reviewed the EHC and lEC discovery potential for leptophilic interactions, which are 
characterized by their negligible coupling to quarks and EW bosons [§]. In paiticular, we have 
discussed in detail the case of an extra gauge boson coupling to E^ — E^, as it also provides a 
good estimate of the reach of these machines. If it only couples appreciably to muons and taus, 
as assumed, it can be only observed as a resonance in the di-muon distribution of four-lepton final 
slates. Then, fhe corresponding Monte Carlo simulafion shows fhaf fhe besf EHC (lEC) limifs can 
be derived using 3p plus missing energy {4p) evenfs. We have found lhal fhe EHC can be sensilive 
up lo Z'^_^ masses near fhe TeV for g' of order 1. The lEC can be only sensilive lo Z'^_^ masses 
below 300 GeV bul if could evenlually provide a direcl measuremenl of fhe veclor boson couplings 
fo tau leptons. Similarly, we have sfudied fhe corresponding bounds for fhe ECC, finding fhaf if is 
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sensitive to gauge boson masses of ~ 2.5 TeV for g' ~ 1. and to gauge couplings almost 5 times 
smaller than at the LHC for lower Z'^_^ masses. 

All these results, however, assume a vanishing kinetic mixing of the new vector boson with 
the SM (hypercharge). Otherwise, the current coupling to the extra vector boson gets a new term 
proportional to the kinetic mixing and to the hypercharge and hence, the new vector boson also 
couples to electrons and to quarks but with couplings suppressed by the kinetic mixing and ratios 
fixed by the SM hypercharges. What implies new contributions to Drell-Yan production mediated 
by the new gauge boson, improving appreciably the LHC and ILC discovery potential due to their 
sensitivity to lepton pair production, pp —)■ 7 ,Z,Z' —ll and e^e~ —)■ 7 ,Z,Z' —)■ ll, respectively. In 
fact, the new contributions are in both cases proportional to the kinetic mixing square. Limits from 
LEP 2 force this mixing to be at most of few per cent for relatively low Z^_.j. masses (~ 300 GeV) 
but leave it essentially unconstrained for gauge boson masses of the order of few TeV. In contrast, 
the LHC data already provide the strongest bounds on this mixing, reaching the per mille level if no 
departure from the SM is observed with the full expected energy and luminosity. This limit could 
be eventually improved by another order of magnitude at the FCC. 

Few comments are in order to conclude. This study enforces the need to give experimental 
results specifying the flavor in the sampling, and to provide limits as model independent as possible 
on the corresponding cross-sections [|^]. Although the kinetic mixing may receive small quantum 
corrections, its size can be large at low energy due to the matching with a more fundamental 
theory at high scale, hence making the phenomenological analysis compulsory. Finally, the SM 
extension with an extra gauge boson coupled to F^ — F^ has been recently invoked to explain small 
discrepancies in Higgs and B decays p^, although the proposed models include extra matter and 
assumptions beyond the phenomenological analysis performed here. Feptophilic interactions have 
been also considered to produced and annihilate dark matter In any case, many approaches 
corner NP in the heavier families, especially when dealing with the lepton sector 
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